Oscillating Bispectra and Galaxy Clustering: A Novel Probe of Inflationary Physics 

with Large-Scale Structure 
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Many models of inflation predict oscillatory features in the bispectrum of primordial fluctuations. 
Since it has been shown that primordial non-Gaussianity can lead to a scale-dependent halo bias, 
we investigate the effect of oscillations in the three-point function on the clustering of dark-matter 
halos. Interestingly, we find that features in the inflaton potential such as oscillations or sharp steps 
get imprinted in the mass dependence of the non-Gaussian halo bias. In this paper, we focus on 
models displaying a sharp feature in the inflaton potential as well as resonant non-Gaussianity. In 
both cases, we find a strong scale dependence for the non-Gaussian halo bias with a slope similar 
to that of the local model. In the resonant case, we find that the non-Gaussian bias oscillates with 
halo mass, a novel feature that is unique to this type of models. In the case of a sharp feature in the 
inflaton potential, we find that the clustering of halos is enhanced at the mass scale corresponding to 
the Fourier mode that exited the horizon when the inflaton was crossing the feature in the potential. 
Both of these are new effects that open the possibility of characterizing the inflationary potential 
with large-scale-structure surveys. We briefly discuss the prospects for detecting these non-Gaussian 
effects. 



PACS numbers: 98.80.-k,98.80. Jk 



I. INTRODUCTION 



Since the publication of the seminal inflation papers 
P, [1[ , a plethora of models have been proposed to ex- 
plain why the Universe underwent a phase of exponential 
expansion at early times. Since most models offer very 
similar basic predictions, distinguishing between these 
models with today's data is not an easy task. One ap- 
proach that has received a lot of attention recently is 
to look for departures from Gaussianity in the primor- 
dial cosmological perturbations [|| . Indeed, while a large 
class of models predicts that the non-Gaussian signature 
should be undctectably small, there also exist a numbcr 
of models for which departures from Gaussianity should 
be relatively large and observable Q. Thus, any detec- 
tion (or absence thereof) of non-Gaussianity in the pri- 
mordial spectrum of perturbations could then rule out a 
large swath of inflation models. 

Non-Gaussian signatures have been traditionally 
looked for in cosmic microwave background (CMB) 
anisotropies [f|. However, it has recently been shown that 
the initial departure from Gaussianity could be amplified 
in the cluste ring of dark-matter halos (see Q for a 
review and [K| for similar effects in another context). 
Indeed, mode coupling in non-Gaussian models induces 
a dependence of the local power spectrum on the long- 
wavelength potential perturbations. This can lead to a 
scale-dependent halo bias on large scales which is observ- 
able in large-scale-structure surveys, since galaxy cluster- 



ing is closely connected to halo clustering on large scales. 
Competitive upper limits on non-Gaussianity have al- 
ready been placed using this method fill ] . 

While the non-Gaussian bias correcti on g oes as k~ 2 in 
the local model, it has been shown in |l2l. Il3j that the 
scale dependence of other types of non-Gaussian mod- 
els can be significantly different. Furthermore, models 
could also differ by how the bias varies with halo mass. 
Therefore, measurements of the biasing of dark-matter 
halos could be used to distinguish among different non- 
Gaussian scenarios. So far, the bispectrum shapes for 
which large-scale-structure predictions have been worked 
out include the local equilateral [lH, orthogonal 

[lH, and folded 17 1 shapes, all of which are scale inde- 
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pendent. However, there are several classes of inflation- 
ary models which predict bispectra that have strongly 
scale-dependent oscillatory features (l7l - [27T | . These mod- 
els can circumvent the tight limit on the bispectrum in 
the squeezed configuration [28| by breaking the slow-roll 
approximation. Since the squeezed triangle configura- 
tion is what determines the scale-dependent halo bias, 
such models potentially leave interesting signatures in 
halo clustering. The oscillatory bispectrum shapes arc 
generally nonfactorizable and are therefore very compu- 
tationally intensive to constrain with CMB data alone 
[2!||3(|. We show here that these models can also be con- 
strained by calculating their impact on halo clustering. 
Moreover, these models leave distinct features in the mass 
dependence of the non-Gaussian halo bias, which allow 
us to distinguish them obscrvationally from the smooth, 
scale-invariant shapes considered thus far. 

In this paper, we calculate the non-Gaussian correction 
to the dark-matter halo bias for two different oscillatory 
bispectra. We focus on models that display a sharp fea- 
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ture in their inflationary potential as well as models that 
have periodic features in the potential. Models with fea- 
tures in the potential have been invoked to explain de- 
viations in the observed CMB power spectrum from the 
smooth prediction. On the other hand, periodic modula- 
tions of the pot ential are motivated by axion-monodromy 
models [24l |25|. While in both cases one obtains oscil- 
latory three-point functions, the physics responsible for 
these modulations is very different. Indeed, the non- 
Gaussianities in the model with a feature are generated 
when the mode exits the horizon while for the resonant 
model, the non-Gaussianities are generated deep inside 
the horizon. As a consequence, we expect the two infla- 
tionary scenarios to make distinct predictions about the 
clustering of dark-matter halos. In particular, we antic- 
ipate that in the feature model, the non-Gaussian effect 
should be the largest around the mass scale that exited 
the horizon while the inflaton was crossing the feature. 
On the other hand, we expect non-Gaussian effects to 
be important for a broad range of scales in the resonant 
model since these were generated by causal physics in- 
side the horizon. Our results support these qualitative 
predictions and most interestingly, they allow us to map 
properties of the inflaton potential to features of galaxy 
clustering. 

The structure of this paper is as follows. Wc begin 
by briefly reviewing halo biasing in the peak-background 
split formalism for general nonlocal quadratic non- 
Gaussianity. We then calculate the scale-dependent cor- 
rection to the halo bias for the two inflation models con- 
sidered here, emphasizing the effect of the new term un- 
veiled in [HI, H3 • We finally discuss our results in light 
of our qualitative predictions and physical expectations 
and conclude with a discussion on how these new effects 
could be detected in large-scale-structure data. 



is usually parametrized by the bispectrum, 

($(k 1 )$(k 2 )$(k 3 )) = (2^) 3 ,5| 5 (k 1 +k 2 + k3) J B 4 ,(fc 1 ,fc 2 ,fc3): 

(2) 

where 3d is the Dirac delta function and k 3 = ki + 
k 2 |. In terms of the Fourier space kernel W(ki, k 2 ), the 
bispectrum amplitude is given by 



B*(fci,fc2,fc 3 ) = 2/nl W(ki,k 2 )P*(A;i)P*(fc 2 ) 



2 perm. 



(3) 



where P$(fc) stands for the power spectrum of $. The 
two permutations not written are the two remaining 
cyclic permutations of ki, k 2 , k 3 . Since the kernel 
W(ki,k 2 ) is only required to be symmetric under the 
exchange of its two vectorial arguments, Eq. ([3]) does 
not uniquely specify W. One possible choice of kernel is 



W(ki,k 2 



1 



B^(k ll k 2 ,k 3 ) 



2/nl P*(fci)P#(fei) + 2 perm. 



(4) 



which has the nice property of being fully symmetric un- 
der the exchange of the three momenta. For the halo 
bias calculation, we are mainly interested in the squeezed 
limit of the kernel where k 2 ,k 3 3> k\. In this limit, the 
bisp ectrum uniquely defines the kernel via the relation 

m 



wcki,k 2 ; 



fe 2 >fci p«t.(fci, fc 2 , fc 3 ) 

^ 4/ N LP*(fc 1 )P$(fc 2 )' 



(5) 



Finally, to compute the dark-matter halo bias at 
late times, we need to consider the processed kernel 
Wo(ki,k 2 ) defined via the transfer function T(k), 



II. NON-GAUSSIAN HALO BIAS IN THE 
PEAK-BACKGROUND SPLIT FORMALISM 

A. Nonlocal Kernel 

Following [fJJ], we consider the case for which the 
Bardeen potential during matter domination $ is a gen- 
eral, nonlocal quadratic function of a Gaussian field </> 

$(x) = 0(x) + / NL J d 3 y J d 3 z^(y,z)^(x + y)0(x + z), 

(1) 

where the kernel W(y, z) is symmetric in its arguments 
and only depends on t/, z as well as y • z. In Fourier 
space, one can think of F7(ki,k 2 ) as a scale-dependent 
coupling between different modes. To conform to stan- 
dard notation, we have pulled out an arbitrary factor of 
/nl from the non-Gaussian kernel. Our results do not de- 
pend on this particular choice as they are only sensitive 
to the product fni,W(y,z). Deviation from Gaussianity 



W (ki,k 2 ) 



Tflki+kaDW^ki.kz) 



T(h)T(k 2 ) 
In the squeezed limit, this reduces to 



W (ki,k a ) 



T(fci)4/ NL P*(fci)P*(fc 2 )" 



(6) 



(7) 



Note that we define our /nl in terms of the Bardeen 
potential at last scattering, conforming to the convention 
usually adopted in CMB analyses. 



B. Halo Bias in Peak-Background Split 

We work in the Lagrangian picture of halo biasing 
where halos are identified as high-density regions in the 
initial linear matter field. As such, we focus here on de- 
riving the Lagrangian halo bias 6/ which relates the halo 
power spectrum to the linear matter power spectrum, 
Ph(k) = bjP(k). The latc-timc linear Eulerian bias rel- 
evant for observations on large scales is simply given by 
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bf = 1 + bj. In the Lagrangian picture, the number den- 
sity of halos per unit logarithmic mass (also called halo 
mass function) is sensitive to the statistics of small-scale 
perturbations. In the Gaussian case, each Fourier mode 
evolves independently and therefore the small-scale mat- 
ter power spectrum P(k a ) (at some initial early time) 
is the same everywhere. However, non-Gaussianity in- 
troduces mode coupling resulting in a dependence of the 
small-scale power spectrum on the local value of long- 
wavelength fluctuations. Non-Gaussian initial conditions 
thus generally rescale the local small-scale variance of the 
density field smoothed over a scale R s , eros, according to 

m 



'Os 



4/NL<Mk)criU£;), 



where the spectral moment a^y is 
d 3 k 



F^(k s )W (k,k s )P(k s 



(8) 



(9) 



and where 4>l (k) is a long- wavelength fluctuation of the 
gravitational potential. Here, Fr s is the Fourier trans- 
form of a spherical tophat with radius R s , P(k B ) is the 
matter power spectrum, and the "hat" denotes quanti- 
ties that contain non-Gaussian contributions. Notice the 
appearance of the non-Gaussian kernel which indicates 
how a mode with wave number k s couples to the long- 
wavelength mode k. Note also that a^ is not positive 
definite, as the sign depends on the shape of the non- 
Gaussian kernel. However, the second term in Eq.[S]is al- 
ways much smaller than the Gaussian variance a$ s (since 
4>l ~ 10~ 5 ), so that a 2 s is always positive. 

Since the halo abundance hh generically depends on 
<7os, this induces a scale-dependent dark-matter halo bias 
of the form [H El 



b!{M,z;k) 



1 dn h 



n h d5 L (k) 
<91nn/j 
d\np din oos dSL{k) 



S L =Q 

d In fi/j d In ctq s 



(10) 



where M stands for the halo mass, z for rcdshift, p is 
the average matter density of the Universe, fih is the 
average number density of halos of mass M and Si(k) 
is a long-wavelength density fluctuation. The halo mass 
M is related to the smoothing scale R s through M — 
(4ir/3)pRs for a spherical tophat window function. In 
the following, we will drop the explicit z dependence. 
The first term in Eq. (fTU|) is the usual Gaussian bias b\ 
while the second term is induced by the non-Gaussian 
initial conditions. This last term can be expressed in a 
compact way when adopting a universal mass function 
prescription, 



M 



d In &o s 



dlnM 



(11) 



where v = 5 c /ao s is the significance, S c m 1.686 is the 
linearly extrapolated collapse threshold, and f(v) is a 



multiplicity function which we do not need to specify 
explicitly. A change in <7q s thus changes halo abundance 
through a change in v as well as a change in the Jacobian 
|91n<xo s /91nM|. The non-Gaussian halo bias correction 
can then be written in terms of the non-Gaussian kernel, 



Ab I (M,k)=2fa h M- 1 (k) 



' w 



(M,k) 



with 



e w (M,k) 



x [b 1 (M)5 c + 2e w {M,k)} 



d\na 2 w {M,k) 



dhxa 2 s {M) 



(12) 



(13) 



where M(k) = 2k 2 g*(z)/(3(l + z)H 2 Q m ). Here g*(z) 
is the potential growth function normalized to unity at 
last scattering. Since it is understood that k s 3> k in 
Eq. ([5]), we see that the bias correction depends on the 
non-Gaussian kernel evaluated in the squeezed limit. We 
note that the term proportional to e^(M, k) had been 
previously neglected in the literature until it was shown 
to be important in [3lj | . As we will see in the next sec- 
tion, this term is crucial for models displaying oscilla- 
tory features in their bispectrum. Examining Eq. (|12l) . 
we observe that the scale dependence of the halo bias is 
determined by the product of Ai~ 1 (k) oc k~ 2 with the 
leading fc-dependent part of the processed non-Gaussian 
kernel evaluated in the squeezed limit. We now turn our 
attention to bispectra showing oscillatory behavior and 
calculate the resulting scale-dependent bias. The numer- 
ical results presented in this paper assume a flat ACDM 
universe with h = 0.72, Q m = 0.28, n s — 0.958 and 
as = 0.8. The pivot scale for the primordial power spec- 
trum amplitude is kept at k„ = 0.002Mpc~ 1 throughout. 



III. OSCILLATORY BISPECTRA AND THEIR 
SCALE-DEPENDENT BIAS 

A. Resonant Non-Gaussianity 

Resonant non-Gaussianity arises when periodic fea- 
tures in the inflationary potential lead to an oscillatory 
coupling between modes, which can trigger a resonance 
for modes oscillating with the same frequency inside the 
horizon (2^, Hi], I23H251 [27j ■ Such features arise, for ex- 
ample, in certain brane inflation models or in axion- 
monodromy inflation. For this class of models, the bis- 
pectrum has the generic form (27j 



*«»=(!) ( 2 ^) 4 /NL S A|^p2 y S m{C u \n{k t /k p )) 
+ i- cos (C u In (k t /k p )) £ | + O (J^j ^ , 



(14) 
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where A<j, is the amplitude of primordial scalar power 
spectrum, k t = k 1 + k 2 + k 3l k p is a pivot scale which 
introduces a phase, and C u is related to the frequency w 
of the periodic features of the inflationary potential by 
C w = uj/Hi. Here, Hi stands for the Hubble parameter 
during inflation. The leading factor of 5/3 comes from 
the conversion between the Bardeen primordial poten- 
tial $ and the gauge invariant curvature perturbation £ 
at late times. Constraints from the matter power spec- 
trum provide an upper bound on the value of f^. For 
the axion-monodromy scenario with a linear zero-order 
potential, this bound reads [27], H(| 

/nT < W- 3 C 5 J 2 , (15) 

where it is assumed that the pivot scale k p = 0.002 
Mpc -1 exits the horizon about 60 e- folds before the end 
of inflation. Other zeroth-order inflaton potentials are 
likely to lead to a somewhat different constraint on 
but we shall use Eq. (Tl3|) as a rough upper limit for this 
type of model. The resonant bispectrum can readily be 
evaluated in the squeezed limit [27( 
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FIG. 1: Non-Gaussian correction to the halo bias for the 
resonant non-Gaussianity model as a function of scale. We 
evaluate the bias for M = 10 13 M Q //i at z = 0. We take 

/nl = W~ 3 C^ 2 and evaluate the Gaussian bias bi using the 
Sheth-Tormen mass function 1331. 



B T , 



fc s »/c 



res o A 2 
NL 



C u In 



2 ha 



(16) 

Note that the leading correction to this expression is sup- 
pressed by a factor of k/k s , thus negligibly contributing 
to the bispectrum in the squeezed limit. Using Eq. (|5|), 
we obtain the leading-order non-Gaussian kernel, 



VF(k,k s ) 



(2tt) 4 ( k 

k 
k 



~k 



2fc, 



2 V fc s 



3-e 



. (17) 



where e = n s — 1 and the ellipsis stands for terms that 
are higher order in k/k s . We immediately see that the 
scale dependence of the non-Gaussian bias is given by 



A6 /jres (fc) oc kT 



(18) 



that is, it is very similar to that of the local model. To 
calculate the amplitude of the bias correction, we first 
need to integrate Eq. © over the small-scale modes to 
obtain the non-Gaussian spectral moment a 2 v (k). The 
k s integral is of the general form 



aw oc / k s dk s T 2 (k s ) 



]l{k s Rs) 

m 



cos (C w In (2k s /k p )). 



(19) 

For large values of the frequency C w , the integrand is 
rapidly oscillating and the resulting amplitude for the 
non-Gaussian bias is expected to be rather small. For 
a small enough value of the frequency (C^ < 100), the 
integral can be done numerically. To evaluate the second 
term in Eq. (|12l) , we first use the chain rule to write it as 



da 2 w (M,k) 
da 2 0s (M) 



da 2 w (M,k) 
d\nM 



dlnM 



(20) 



The derivatives on the right can be calculated numeri- 
cally. For all numerical computations, we use the com- 
plete expression for the bispectrum, Eq. (fl4|). In Fig. [U 
we show the scale dependence of the non-Gaussian halo 
bias correction for three values of C w evaluated for a 
halo mass of 10 13 M Q /h at z = 0. We see that the non- 
Gaussian bias correction is small except for the largest 
scales where the scale dependence of A6/ ires (fc) oc k~ 2 ~ e 
becomes important. Interestingly, the non-Gaussian bias 
for these resonant models is completely dominated by the 
term proportional to ew{M, k) in Eq. (|12p . which was re- 
cently unveiled in [3ll |32| . To understand why this new 
term is crucial for our analysis, we plot in Fig. [2] the non- 
Gaussian spectral moment a^y as a function of halo mass 
for a fixed comoving scale. We see that a\y strongly os- 
cillates with halo mass, leading to a large contribution to 
da^/d In M, especially toward small masses. This high- 
lights the importance of the newly discovered term for 
accurately predicting the non-Gaussian halo bias. We 
will discuss the relevance of this result for observations 
in Sec. Hvl 



An interesting feature of resonant non-Gaussianity 
models is that they predict a modulation of the halo bias 
with changing halo mass. In Fig. [31 we show the non- 
Gaussian halo bias as a function of halo mass evaluated 
at a scale k = 10~ 3 /i Mpc -1 . Again, the bias is domi- 
nated by the second term of Eq. |12]) for M < W 15 M Q /h. 
We observe that the amplitude of the non-Gaussian bias 
decreases with increasing C u very rapidly and therefore 
this effect is likely to be unobscrvable unless is small. 
As expected, the non-Gaussian features of the halo bias 
show coherent modulations over a wide range of mass 
scales, an artifact of non-Gaussianitics being produced by 
causal physics deep inside the horizon for these models. 
To contrast the resonant model with the more traditional 
local model of non-Gaussianity, we also display the halo 



5 



— c u . 


= 10 


— Co, 


= 20 




FIG. 2: Non-Gaussian spectral moment a%/ for the resonant 
model as a function of halo mass. We evaluate this spectral 
moment for k — 10~ 3 h Mpc -1 at z = 0. 
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FIG. 3: Non-Gaussian correction to the halo bias for the res- 
onant non-Gaussianity model as a function of halo mass. We 
evaluate the bias for k = 10~ 3 h Mpc -1 at z = 0. We take 
/nl — 10~ z Ct/ 2 and evaluate the Gaussian bias b\ using the 
Sheth-Tormen mass function. For comparison, we also show 
the bias for local non-Gaussianity with /Jj°L al = 2. 



the authors of [2Cj suggested an approximate form: 

^ff^^k (21) 
k t Ak f 



k-k 2 



sin 



sin 



k 2 
K f 



Here, k t = k\ + k2 + k% and kf ± Akf are the Fourier 
modes that exit the horizon while the inflaton is crossing 
the feature (the sharper the feature in the inflaton poten- 
tial, the larger Akf becomes). We fixed the phases such 
that -Bfeat - > as k t — >• which is physically motivated 
since modes that exit the horizon long before the infla- 
ton encounters the feature should not show significant 
non-Gaussianities. We choose the overall sign such that 
the non-Gaussian bias is positive for the scale exiting the 
horizon when the inflaton crosses the feature. In prac- 
tice, this sign should be fixed by comparison to numerical 
simulations. In the squeezed limit, the bispectrum reads 



(27T) 4 /£f 



A 2 
kik 2 



x sm 



1k s 



2k, Ak 



sm 



(22) 



For a narrow feature, we generally expect Akf/kf <C 1 
and thus the last sinusoidal factor in Eq. (|22j) can be con- 
sidered as an envelope function for the first rapidly os- 
cillating sine factor. From Eq. (|22|) . we see that at fixed 
k s , non-Gaussianity becomes more important for modes 
k smaller than kf (remember that k is the scale at which 
clustering of halos is measured). The non-Gaussian ker- 
nel can then readily be obtained as 




x sm 



(23) 



bias for a local model with /^L al = 2. At moderate halo 
masses (M ~ 10 12 — 10 14 M Q //i), the effect of resonant 
non-Gaussianity is comparable to that of a local model 
with /^L al of order unity, for the values of C w chosen. 



B. Features in the Inflaton Potential 

The presence of a sharp feature in the inflaton poten- 
tial can induce large primordial non-Gaussianities [l8l - 
|22T | . Indeed, modes that exit the horizon while the infla- 
ton is crossing the feature get a boost in their three-point 
signal. Here, we shall focus on the case of a step in the 
inflaton potential, but our analysis could also be applied 
to the case of a bump in the potential. The exact form 
of the bispectrum can only be obtained numerically but 



Note that W approaches zero for k s <C kf and for 
k s 3> kf. This physically makes sense since modes with 
k s 3> kf are oscillating deep inside the horizon when 
the inflaton crosses the feature and we expect their non- 
Gaussianities to roughly cancel out. On the other hand, 
modes with k s <Si kf are outside the horizon when non- 
Gaussianities are generated and we thus expect their con- 
tribution to the kernel to be small. We immediately see 
that the scale dependence of the halo bias is given by 



A6/,feat(fc) oc k" 



(24) 



which, at first look, is similar to the resonant model. 
However, as we will see below, the two models predict 
very different behaviors for how the amplitude of Abi 
varies with halo mass. From Eq. (|23|) . it is straightfor- 
ward to compute numerically the non-Gaussian spectral 
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FIG. 4: Absolute value of the non-Gaussian spectral moment 
o^y for the feature model as a function of halo mass. We 
evaluate this spectral moment for k = 10 _3 /i Mpc -1 at z = 
and use Akf/kf = 0.01. 

moment crj^. In the feature model, the integral over small 
scale modes has the general form 

2 f u „ rtn. , jl(k s Rs) fk f \ 
a w oc / k s dk s T (k s ) — I I 

^HT} (25) 

We see that for a given value of kf, there will always be a 
scale i? s for which constructive interference between the 
first sine factor and the Bcsscl function happens (again, 
the second sine factor is considered to be slowly varying). 
Since the Bessel function peaks around k s ~ R ~ , we 
naively expect \&yy \ t° nave a maximum near i? s = Rf ~ 
0(l)/kf or equivalents, near M/ - e>(l)(47r/3)p£^ 3 . 
For i? s <C i?/, the first sine factor is rapidly oscillating 
near the peak of the Bessel function and thus the only 
nonvanishing contribution comes from the low k s tail. 
In this limit, the integrand becomes independent of R s 
since jf(k s R s ) oc k 2 s R 2 s for k s <C R^ 1 . We thus expect 
a w to asymptote to a constant for small halo masses. On 
the other hand, for R s 3> Rf, the integrand approaches 
zero and we therefore expect a w to vanish for large halo 
masses. 

Our numerical calculations support these qualitative 
conclusions. Indeed, Fig. @] shows the absolute value 
of the non-Gaussian spectral moment for two values of 
kf. The enhancement around Mf ~ 1.5(4n/3)pkJ 3 is 
clearly visible in both cases and a w (M) rapidly vanishes 
for M > Mf. The latter point is expected on physical 
grounds since halos with M > Mf correspond to scales 
(in the initial density field) that were outside the horizon 
when the non-Gaussianities were generated. Therefore, 
the non-Gaussian correction to the variance of the den- 
sity field smoothed on these scales must be somewhat 
suppressed. 

We can now use Eq. (jT2J) to compute the scale- 




M [M Q /h] 

FIG. 5: The two contributions to the non-Gaussian halo bias 
correction as a function of mass for the model with a feature 
at k f = 0.5hMpc _1 . We evaluate a\ r for k = 10~ 3 hMpc _1 
at z — and use Akf/kf = 0.01. The Gaussian bias bi is 
derived using the Sheth-Tormen mass function. 



dependent non-Gaussian correction to the halo bias. 
From the functional form of a w , we expect the term 
proportional to ew to dominate around the feature at 
71/ = Mf since d In a w jd In M is largest there. To verify 
this, it is instructive to consider the two distinct contribu- 
tions to the bias as a function of halo mass. In Fig. [3J we 
show both the contribution proportional to the Gaussian 
bias b\ as well as the recently unveiled contribution pro- 
portional to ew for a feature at kf = 0.5 hMpc^ 1 . Wc 
observe that the ew term clearly displays a feature at 
Mf - 1.5(47r/3)pfcJ 3 ~ 4 x 10 12 M Q /h and that it dom- 
inates the overall bias for halo masses M < 10 14 M Q //i. 
This once again highlights the importance of the ew term 
for models with strongly scale-dependent bispectra (3^ |. 

In Fig. [6l we show the amplitude of the non-Gaussian 
bias correction as a function of halo mass for three dif- 
ferent values of kf. Note that while we assumed a 
value of /nl* = ^0 here, still within the allowed range 
of power spectrum contraints [2fj| . this parameter is in 
reality fixed for a given inflation potential by the nu- 
merical calculation. For comparison, we also plot the 
halo bias correction for local quadratic non-Gaussianity, 
Afr/jocai = 2f^ aX bi8 c M~ 1 (k). We immediately see that 
the bias correction for models with a feature in the po- 
tential displays a large enhancement around M = Mf 
when compared to the monotonic and featureless bias 
of a local-type model. Observations of this tell-tale sig- 
nature in large-scale-structure data provide us with an 
exciting new window to probe microscopic inflationary 
physics with observations of galaxy clustering on the 
largest scales. 

For halo masses M ^> Mf, the bias becomes domi- 
nated by the first term of Eq. (fT2|) since b\ is large for 
very massive halos. In this limit, Abj is very sensitive to 
the small k behavior of the bispectrum which may not be 
accurately captured by our ansatz Eq. (|22p . Thus, a com- 
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FIG. 6: Mass dependence of the non-Gaussian correction to 
the halo bias for the model with a feature in the inflaton 
potential. We evaluate the bias for k — 10 _3 /i Mpc -1 at 
z = and use Ak f /k f = 0.01. We take = 10 and 

evaluate the Gaussian bias b\ using the Sheth-Tormen mass 
function. For comparison, we also show the bias for the local 
model of non-Gaussianity with /JjL al =-10. 
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FIG. 7: Scale dependence of the non-Gaussian correction to 
the halo bias for the model with a feature in the inflaton 
potential. We evaluate the bias for M — 1.5(47r/3)pfcJ 3 at 
z = and use Akf/kf = 0.01. We take /£ e £* = 10 and 
evaluate the Gaussian bias 61 using the Sheth-Tormen mass 
function. For comparison, we also show the non-Gaussian bias 
correction from the local model evaluated at M = 10 13 Mq /h 
and /J,°£ al = 100. 



plete numerical computation of the bispectrum is likely 
to be required to accurately predict the large-mass limit 
of the halo bias. We leave this for future work. However, 
since halos above M > 10 15 M Q /ft, are very rare, espe- 
cially at higher redshifts, we do not expect the observa- 
tional constraints to be dominated by this mass range. 

In Fig. we show the absolute value of the non- 
Gaussian bias correction as a function of scale for three 
values of kf. We evaluate the amplitude of the bias at 
AI = Mf, that is, at the peak of the feature in ew- At 
this mass scale, we see that the halo bias can reach a very 
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wide amplitude on large cosmological scales. The scale 
dependence proportional to k~ 2 ~ e derived in Eq. is 
readily visible. For comparison, we also plot the halo bias 
correction for the local-model of non-Gaussianity which 
displays a similar scale dependence but a much smaller 
amplitude even for /^L al = 100. 

We note in passing that varying Akf/kf corresponds 
to changing the overall scale of the bispectrum. This 
can readily be seen from Eq. (|22[) where Akf/kf only 
appears in the modulating envelope. For Akf/kf <C 1 
and 2k s < kf (i.e. for the modes that contribute most 
to ), one can Taylor expand the second sine factor to 
obtain, after simplification, -Bf ca t oc Akf/kf. Therefore, 
the overall amplitude of the bias correction is determined 
by the product flf^(Akf /kf). This scaling agrees with 
the result of [19j . where it was shown that the overall 
amplitude of the bispectrum is inversely proportional to 
the width of the step in the potential. To see this, we 
note that the sharper is the step, the more kinetic en- 
ergy is acquired by the inflaton and by consequences, the 
longer slow roll is violated. Therefore, we expect that the 
sharper is the step (corresponding to larger bispectrum 
amplitude), the larger the band of Fourier modes affected 
(Akf) will be, hence the above result. 

In summary, we have shown that the presence of a 
feature in the inflaton potential leads to a correspond- 
ing feature in the mass dependence of the non-Gaussian 
halo bias. Ultimately, this is a consequence of non- 
Gaussianity being generated at a specific scale during 
inflation in these models. Finally, we reiterate that the 
numerical results presented in this section were computed 
using the analytical expression for the bispectrum given 
in Eq. 1211 It is important to keep in mind that this ex- 
pression is approximate. However, it most likely captures 
the important physics. As such, we expect our conclu- 
sions to be robust to the inclusion of a more accurate 
bispectrum. 



IV. DISCUSSION 

We have analyzed the non-Gaussian halo bias re- 
sulting from two inflation models displaying oscillatory 
bispectra. Even though the two models predict the same 
scale dependence as local quadratic non-Gaussianity, we 
find that they make very different predictions concerning 
how the amplitude of the bias varies with halo mass. 
Indeed, while the resonant model predicts an oscillatory 
amplitude as a function of halo mass, models with a 
feature in the potential predict an enhancement of the 
bias for halos with mass that corresponds to the scale 
that exited the horizon at the time when the inflaton 
was crossing the feature in the potential. Ultimately, 
these very different outcomes can be traced back to 
the distinct physics that is responsible for generating 
non-Gaussianities in the first place. 

For the resonant model, non-Gaussianity is generated 
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well inside the horizon when the modes are rapidly 
oscillating. As explained in J2(j, oscillations in the 
inflaton potential lead to an oscillatory coupling between 
different Fourier modes. As the physical frequency of 
each mode k/a(t) decreases, there will be a time when 
k/a(t res ) ~ u! and the oscillating mode can resonate with 
the coupling and generate a departure from Gaussianity. 
As a large number of modes eventually passes through 
the resonance, we naturally expect the non-Gaussian 
effects to be present on a broad range of scales, as can be 
seen in Fig. [3] A crucial consequence of this subhorizon 
generation mechanism is that the resulting modulation 
of the halo bias is in phase over a broad range of masses. 
This is a tell-tale signature that could be looked for in 
large-scale-structure data and used to put constraints 
on these resonant-type inflation models. 

For models with a feature in the inflaton potential, 
non-Gaussianity is generated during slow-roll violation 
associated with the inflaton suddenly accelerating as it 
crosses the step or the bump in the potential. Conse- 
quently, different Fourier modes within a limited range 
of scales are coupled, hence generating a nonvanishing 
three-point function. However, modes deep inside the 
horizon are rapidly oscillating and we thus expect their 
non-Gaussian signature to average out to zero. On the 
other hand, modes that exit the horizon as slow roll is vio- 
lated are frozen-in before causal physics could erase their 
correlation with other Fourier modes. We thus expect the 
bispectrum to be significant when at least one side of the 
triangle has k ~ kf (and no side with k 3> kf). As a con- 
sequence, modes that exit the horizon during slow-roll vi- 
olation get an enhanced coupling to the long-wavelength 
perturbations resulting in an amplified clustering of ha- 
los at the corresponding mass scale. Conversely, modes 
that are superhorizon when slow-roll is violated become 
correlated with modes that have k ~ kf. This induces 
a rescaling of the variance of the density field according 
to Eq. ([SI) which results in a nonvanishing halo bias at 
these mass scales. As mentioned earlier, this rescaling 
of the variance is very sensitive to the small-k limit of 
the bispectrum and a complete numerical computation 
will be required to accurately predict the halo bias for 
M 3> Mr. Nevertheless, since very massive halos are 
rare and restricted to low rcdshifts, it is unlikely that 
observational constraints will depend sensitively on the 
high-mass tail. 

For both classes of models, we find that the term com- 
ing from the Jacobian dlnaos/dlnM (see [3l[ for details) 
is an important and often dominant contribution to the 
non-Gaussian halo bias. This term comes about because 
the non-Gaussian mode-coupling induces a modulation 
with (f>L of the significance interval dhiv = dlneros 
that corresponds to a fixed logarithmic mass interval 
dlnM. This term is strictly present because we have 
assumed throughout that halos are selected by mass, 
which is appropriate e.g. when comparing to A-body 
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FIG. 8: Comparison between Abi(M) obtained with a spher- 
ical tophat filter and with a Gaussian filter. We evaluate the 
bias for k = 10~ 3 /i Mpc -1 at z = and use Ak//k f = 0.01 
and k f = 0.5ft Mpc -1 . We take = 10. 



simulations. In practice however, galaxies are selected by 
more complex criteria which are only indirectly related 
to the host halo mass. If we divide the total galaxy 
sample into different subsamples (e.g., by luminosity, 
color, light profile, ...), then the mass dependence 
shown in Figs. [3] and [5] will be observable as long as 
the scatter in the mass-observable relation is not much 
larger than the width of the features. Fortunately, the 
latter typically corresponds to a factor of 2 or more 
in mass, which should make these features detectable 
for a wide range of large-scale structure tracers. Note 
that the precise shape of Abi(M) depends on the filter 
chosen for the small-scale density field. In principle, one 
could use a filter matched to the Lagrangian profiles 
of dark-matter halos (34[. We have tried replacing the 
top-hat filter with a Gaussian, and found only relatively 
minor differences, at the 7% level for Abi(M) (see Fig. 
B). 



V. CONCLUSION 

We have shown that measurements of galaxy cluster- 
ing could potentially be used as a probe of features in the 
inflationary potential. By computing the non-Gaussian 
correction to the halo bias, we revealed that features in 
the inflationary potential such as oscillations, bumps, or 
steps get imprinted onto the clustering properties of dark- 
matter halos. While we have restricted ourselves to two 
generic models for which approximate forms of the os- 
cillatory bispectrum are known, we expect this effect to 
be robust to the inclusion of more detailed bispectra. 
We note that this probe of primordial non-Gaussianity 
is complementary to CMB constraints as it probes very 
small scales where the microwave background becomes 
foreground dominated. On intermediate scales, the two 
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approaches could be used in conjunction to cross corre- 
late a possible feature in CMB data with a corresponding 
attribute in the clustering of dark-matter halos. 

While showing the same scale dependence oc fc -2 , the 
predictions of the models considered here are strikingly 
different from the usually considered local model. In 
particular, they show significantly stronger effects for 
moderate halo masses (10 12 — 10 14 Mq/K) than the lo- 
cal model as compared to the effect at the high-mass end 
(> 10 14 M Q /h). Thus, focusing on the most massive, 
highly biased halos might not in general be the best way 
to design or optimize surveys for the search for primordial 
non-Gaussianity. 

The non-Gaussian models discussed here also make 
other predictions which are potentially observable with 
large-scale structure. While the study of these effects is 
beyond the scope of this paper, it would be interesting to 
correlate the non-Gaussian halo bias with features in the 
matter power spectrum and in the mass function of dark- 
matter halos. The former could in principle be probed 
by weak lensing observations, and the latter through the 



abundance of galaxy clusters. The bispectrum of galaxies 
would also be a precise, albeit more complex and compu- 
tationally expensive, approach to testing these inflation- 
ary models. The key advantage of the scale-dependent 
bias is, however, that it is a unique signature of pri- 
mordial non-Gaussianity which is not easily mimicked by 
other effects. We thus anticipate this observable to be a 
robust probe of features in the inflaton potential. 
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